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Abstract We show that all the time-dependent statistical properties of the rightmost points of a branching 
Brownian motion can be extracted from the traveling wave solutions of the Fisher-KPP equation. We 
show that the distribution of all the distances between the rightmost points has a long time limit which 
can be understood as the delay of the Fisher-KPP traveling waves when the initial condition is modified. 
The limiting measure exhibits the surprising property of superposability: the statistical properties of the 
distances between the rightmost points of the union of two realizations of the branching Brownian motion 
shifted by arbitrary amounts are the same as those of a single realization. We discuss the extension of our 
results to more general branching random walks. 
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1 Introduction 

A branching random walk is a collection of points which, starting from a single point, diffuse and branch 
independently of the time, of their positions or of the other points, as in figure [T] 




Position Position 



Fig. 1 Two examples of branching random walks. Left: a branching random walk with discrete time where each point splits into 
two points at each time step. Right: a continuous version called branching Brownian motion where points diffuse as in a Brownian 
motion and branch with a constant rate. 
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Branching random walks appear in many contexts ranging from Mathematics [25 8 22] to Biology 
B16ll20|[T7l . They can for example be used to describe how a growing population invades a new envi- 
ronment. In the one dimensional case, see figure [T] there is, at a given time /, a rightmost individual at 
position X\(t), a second rightmost at Xi{t) and so on. (Note that the rightmost X\{t') at a time t' > t is 
not necessarily a descendant of the rightmost X\ (f ) at time t .) The expected position m t = (X\ (f )) of the 
rightmost individual as well as the probability distribution of its position X\ (t) around m t are well under- 
stood 1 8 25 1 ; the goal of the present paper is to describe the statistical properties of the positions of all the 
rightmost points in the system, in particular the distribution of the distances between the two rightmost 
points, the average density of points at some fixed distance from the rightmost X\ (t ), etc. 

One motivation for studying these distances is that the problem belongs to the broader context of 
extreme value statistics [33 5, 11 24 32 30, 10, 18]: Trying to understand the statistical properties of the 
rightmost points in a random set of points on the line is a problem common to the studies of the largest 
eigenvalues of random matrices [33], of the extrema of random signals [ 12, 13,29,4, 19,6), or of the low 
lying states of some disordered systems such as spin glasses [26,30.7. 1 ,2|. In fact, the points generated 
after some time t by a branching random walk can be viewed as the energies of the configurations of 
a directed polymer in a random medium 11411271 . and the distances between the rightmost points as the 
gaps between the low lying energy states. 

The most studied example of branching random walk is the branching Brownian motion: one starts 
with a single point at the origin which performs a Brownian motion and branches at a given fixed rate 
(right part of figure [TJ. Whenever a branching event occurs, the point is replaced by two new points 
which evolve themselves as two independent branching Brownian motions. While the number of points 
generated after some time t grows exponentially with time, the expected position m, of the rightmost 
point increases only linearly with time [25,8]. In one dimension, Mc Kean ll25l and Bramson [8] have 
shown that the probability distribution of the rightmost point is given by the traveling wave solution of 
the Fisher-KPP equation, with a step initial condition. Here we will see that all the statistical properties of 
the rightmost points can be understood in terms of solutions to the Fisher-KPP equation with appropriate 
initial conditions [9[. We will also show that the distribution of the distances between these rightmost 
points has a long time limit which exhibits the striking property of superposability: the distances between 
the rightmost points of the union of two realizations of the branching Brownian motion have the same 
statistics as those of a single realization. 

This paper is organized as follows: in section[2]we introduce some generating functions useful to study 
random sets of points on the line and show how one can use them to obtain all the properties of these 
random sets. In section[3]we show that, for the branching Brownian motion, all these generating functions 
are solutions of the Fisher-KPP equation. We also show that the distribution of all the rightmost points 
as seen from m t or, alternatively, as seen from Xi (t ), has a long time limit which can be computed as the 
delay of Fisher-KPP traveling waves. This distribution has the property of superposability. In section [4] 
we present results, mostly numerical, on some specific aspects of the limiting distribution of points in 
the branching Brownian motion, namely the distribution of the distance between the two rightmost points 
and the average density seen from the rightmost point. In section [5] we explain how the results on the 
branching Brownian motion can be extended to more general branching random walks. Finally, we study 
in section[6]the distribution of all the rightmost points in a specific frame which depends on the realization 
and which was introduced by Lalley and Sellke E2ll . 

2 Statistics of point measures on the line 

In this section, we introduce some useful quantities (generating functions) to characterize random sets of 
points on the line such that the number n(x) defined as 



n(x) = (the number of points on the right of position x) 



(1) 
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is finite and vanishes for x large enough. 



2.1 The generating functions 

The first generating function one can define is 

= (A"M>. (2) 

From the knowledge of this function, one can extract the probability distribution function pi(x) of the 
position x of the z'-th rightmost point. Indeed, by definition |2]) of i//^, 

= (3) 

;>o 

where Qj(x) is the probability that there are exactly i points on the right of x. One can notice that Qo(x) + 

Gi to H h Qi- 1 (x) is the probability to have less than i points on the right of x. The generating function 

of these sums is, from ([3]), 

j^jWx to = goto A + [Goto + Qi M]A 2 + [Qoto + Qi to + 62 to] A 3 + • • • . (4) 

But Goto +Gi to H 1~ C(— i (x) is also the probability that the z'-th rightmost point, if it exists, is on the 

left of x. Therefore, 

A - 

to = L /><■ WA'', (5) 

1 A 

where p;to dx is the probability that the z'-th rightmost point exists and is in the interval [x,x + dx]. (Note 
that / pi(x)dx < 1 is the probability that there are at least z points on the line.) 

The knowledge of y/^ (x) gives in particular the average distances between the points: from d5l, one 
can see that 

J dxxdxW to = (1 - A) [(Xi) + A (X 2 ) +X 2 (X 3 ) + ■■■], 

= (Xi) — A [(Xi) (X 2 )]-X 2 [(X 2 ) (X,)} - ••• , 

where (Xf) = J xpi{x) dx is the average position of the z'-th point (with the convention that Xj = if there 
are less than z points in the system). Therefore 



/ dxx[d x \jf (x) - d x \if k {x)} = ^(<i/., +1 )A'', 
J i>i 



(6) 



where (rfj = (Xj) — is the average distance between the z'-th and the (z'+ l)-th point. 

To obtain the correlations between the positions of pairs of points, one can start, for y < x, from the 
generating function 

vM-*to = ( A " ( V Cy) )- (7) 

The coefficient in front of A'/i 7 in the expansion of y/)^ in powers of A and fi is the probability that there 
are exactly i points on the right of x and j points on the right of y. As in Q, the coefficient of A'/i 7 in the 
expansion of A/(l — A)xjtt/(1— jti)x (x,y) is the probability that there are less than i points on the 
right of x and less than j points on the right of y, which is also the probability that the z'-th rightmost point 
(if it exists) is on the left of x and the j-th rightmost point (if it exists) is on the left of y. Thus, for y < x, 

T^y jzr^^Ykn ( x >y) = E Pij{x,yWli\ (8) 

1 A 1 jl 
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where Pij(x,y) dxdy is the probability that both the 2-th and j'-th rightmost points exist and lie respectively 
in the intervals [x,x + dx] and [y,y + dy]. 

One can generalize (|2|7|i by defining, forxo > X\ > ■■■ > Xk, the generating functions 



V^Ao h {xo,...,x k ) = {x^~X [Xk] ) 



(9) 



of the numbers n(xo), ■ ■ . ,n(xk) of points on the right of positions xo, ■ ■ ■ ,Xk, and get as in (5 8 1 all the 
higher correlation functions. In that way, all the statistical properties of the measure can be derived from 
the knowledge of the generating functions d9l). 



2.2 The measure seen from the rightmost point 



In the following we will often try to characterize the random set of points as seen from the rightmost 
point (i.e. in the frame where the rightmost point is at the origin). To do so, let us define the generating 
functions of the numbers m(z) of points at the right of z in the frame of the rightmost point. (Note that if 
Xi is the position of the rightmost, then m(z) — n{X\ +z) and one has m(z) > 1 for z < and m(z) = 
forz>0.) 



Xi u ...,t, k {zi,...,Zk) = ( 



(10) 



(As in Q, we assume z\ > zi > ■ ■ ■ > Zk-) These generating functions, as in section 2.1 allow one to 
calculate all the statistical properties of the measure in the frame of the rightmost point (in particular the 
distribution of the relative distances between the points). They can be determined from the knowledge of 
the generating functions W}Lo,...,X k ( x o, ■ ■ ■ i x k) defined in (|9j by 



XXj ,...,x k (zi, • • • ,Zk) = J dx<9 A0 y/Q.A! x k (x,x+zi,...,x+Zk)- 



(11) 



In section [4] we will calculate the density of probability Pn(a) that the two rightmost points are 
separated by a distance a (and that there are at least two points on the line) and the average density p(a) 
at a distance a from the rightmost point. From ( 10 1 one can see that 



Pn(a) = -dadpXiii-a) = ~d a J dx d^d XQ \j/ Qll (x 
Then using that oL (x,x — a) = (d x + d a ) i//b^ (x, x — a), one gets 

Pn{a) = -dg J Axd^Q^x.x-a) 



(12) 



(13) 



By a similar calculation one can show that the average density p(a) of points at distance a from the 
rightmost point is 



p(a) = dl J Axd^o^x.x-a) 



(14) 



2.3 Examples 



We now describe a few examples of such measures. 
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2.3.1 A Poisson process with an arbitrary density r(x) 



Our first example is a Poisson process on the line with a density r(x). We assume that r(x) decays fast 
enough to the right so that a rightmost point exists, and that / r(x) dx — °° so that there are infinitely many 
points on the line. 

By definition of a Poisson process, each infinitesimal interval [x,x + dx] is occupied by a point with 
probability r (x) dx and empty with probability 1 — r(x) dx, and the occupation numbers of disjoint intervals 
are uncorrelated. The probability Qi{x) that there are exactly i points on the right of x is given by 



&(*) 



R(x)'e 



i„-Rx) 



where 



r(z) 



dz. 



From this, we obtain VOL W from (2 3 1 and V / /L|i( x ;.y) from Q in the Poisson process: 
Wx {x) = e-^-^\ WXtl {x,y)=e-^-V R ^-^\ 
Using |6|, the generating function of the average (djj+i) between the z'-th and (/+ l)-th points is 

dx 



e -{l-X)R(x) _ g -R(x) 



(15) 



(16) 



(17) 



The probability distribution function Pn(a) that the distance d\ t % is equal to a and the average density 
p(a) seen at a distance a from the rightmost point are given by 



Pn(a) 



dxr(x + a)r(x). 



,-R(x) 



p(a)= / dxr{x — a)r(x)e 



-R(x) 



(18) 



These expressions can be understood directly from the definition of the Poisson process or, with a little 
more algebra, from ( 13 14 1. One can notice that Pn(a) and p(a) are given by the same expression with a 
replaced by —a and are therefore analytic continuations of each other whenever r{x) is analytic. 



2.3.2 A Poisson process with an exponential density e 



In the special case where the density of the Poisson process is an exponential r(x) = exp(— ax), one can 
simply replace R(x) in the previous expressions by exp(— ax) /a. This gives 



yi(x) =exp 



-{l-X)- 



a 



= exp 



-ax -ay 

-M(l-A) 

a a 



so that from (|6]) 

and thus [9] 
One also has from (fT8 



ln(l-A) 



a 



(di,i+i) 
P n (a) = ae- aa , 



1 

ai' 



p(a) = ae° 



(19) 



(20) 

(21) 
(22) 
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2.3.3 Decorated measures 



Start with a collection of points w,-, distributed according to some measure Vi and, independently for each 
point Ui, replace it by a realization of another measure V2 shifted by We say that the points m, are 
decorated by the measure V2 and call the resulting measure as Vi decorated by V2. 

We assume that Vi and V2 are such that the decorated measure has a rightmost point. If the functions 
VA,M> WXii{ x ^y)' ■ ■ ■ f° r tne measure V2 are known, the decorated measure is characterized by functions 
^n(x,y),-- - given by 



Y x (x) = (Hy x (x-u i )) , 



(23) 



where the average is over all realizations {m,} of the measure Vi. For instance, if V\ is a Poisson process 
of density r(u), then 



f^(x)=]~~[ 1 — r(u)du + r(u) Xj/^ (x — u)du 



exp 



\\jfl [x— u) — l] r(u) du 



^x^y) = ex P 



[vx/i ( x - u ,y - u ) - l ] K«) du 



(24) 



2.3.4 Ruelle cascades 



For a decorated measure where the decoration V2 is a Poisson process of density e~ ax , the average over 
the m,'s in ( |23"j ) leads in general to complicated expressions for Wx(x) or f^ fi (x,;y). The expressions for 
Pn(a) and p(a) are however the same as in ( |22| ) for the pure Poisson process of density e~ ax . In fact, 
all the statistical properties of the distances between the rightmost points are the same as those in the 
exponential Poisson process. 

This can be understood from the following reason: decorating the points u\, ... ... by independent 

realizations of a Poisson process of density e~ ax is equivalent to drawing a single realization of a Poisson 
process of density Y<k e ~ ~ = e~ ax Y.k e<1Uk ^ which is just the same as one realization of a Poisson 
process of density e~ ax shifted by the random variable ln(^ jt e a " A )/a. 

The same argument applies to Ruelle cascades, which can be defined as follows [29,6,7,4]: take an 
increasing sequence of positive numbers OL\ < 052 < • • • an d start with a Poisson process of density e~" lX . 
At each step k > 1, each point in the system is decorated by a Poisson process of density e~ a * x . At step k, 
the measure of points in the system is simply, from the previous argument, a Poisson process of density 
g-a*-* globally shifted by a random variable which depends on the positions of the points at step k— 1. 
Therefore, the statistics of the distances of the rightmost points is the same as for the Poisson process of 
density e~°" kX . 



3 The branching Brownian motion and Fisher- KPP fronts 



3.1 The Fisher-KPP equation 

We are now going to see how the generating functions ( 2|7|9 1 can be determined when the random set of 
points on the line are the points generated at time t by a branching Brownian motion. 

To define the branching Brownian motion we start at time t = with a single point at the origin. 
This point diffuses and branches, and its offspring do the same. After some time f, a realization of the 
process consists of a finite number of points located at positions Xj{t) for i — 1,2,3, .. . Then, during the 
next time interval df <C 1, each point, independently of what the others do, moves a random distance 
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Xi(t + dt) -Xj(t) = Tj,(f)V2df with (rj,(f)) = and (r/,(f) 2 ) = 1, and, with probability df, is replaced by 
two new points located at the same position Xi(t). 

For any function (j) one can define the generating function H^(x,t) by 

Ht(x,t) = (jl ^[x-Xiit)]^, (25) 

where the Xi(t) for i = l,...,N t are the positions of the N t points of the branching Brownian motion at 
time f and (•) denotes an average over all the possible realizations. 

By analyzing what happens during the very first time interval df, one can see that the evolution of 
H$(x,t) satisfies 



H${x,t + &t) = (1 -At)(H^(x-r\y/2&,t)j + AtH^{x,t) 1 . (26) 

The first term in the right hand side represents the motion of the initial point during the first time interval 
df and the second term represents the branching event which occurs with probability df during this first 
time interval. Taking df to zero, one gets 

d t H^=d^H^+Hl-H^, (27) 

which is the Fisher- KPP equation [ 16,21 ,25 1. (The Fisher- KPP equation is often written as d t h = d 2 h + 
h — h 2 , but this is identical to ( |27] i by the change of variable h = \ — H§.) Because there is a single point 
at the origin at time f = 0, the initial condition is simply, from (|25]>, 



H^,(x,0)=(j>{x). (28) 

The generating function (|9]l at time f 

VAo h (xo,...,x k ) = (4 xo) ...l<^) (29) 

can be written, for > z\ > ■ ■ ■ > Zk, as 

W^,...A k (x,x + z h ...,x + z k ) = (]J<l>[x-Xi(t)] \ =H+{x,t), (30) 

where the function (j) (x) is given by 

Hx) = X^ e{x) X 1 l - e{x+zl) ---X l k - e{x+ ^\ (31) 
and where 9 (x) is the Heaviside step function defined by 

Bix)J l for ^°' (32) 
v ' [0 forx<0. 

See figure[2]for the general shape of pT) . 

With the choice (31 1 of (j), the generating function Q and, therefore, all the properties of the point 



measure in the branching Brownian motion at time f can be obtained as solutions of the Fisher-KPP 
equation with the initial condition ( pHj ). 

In the special case k = and = of pTj ), i.e. for the initial condition (j)(x) = 9(x), one gets 

Hg(x,t) = Proba(There are no point at time f on the right of x), (33) 
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<j>{x) 



, ^2 


A 2 


1 




-z x -z 2 



Fig. 2 The function {TT) for k = 2. 



and one recovers the well-known fact II251 I8I that the solution He(x,t) of the Fisher- KPP equation with a 
step initial condition is the cumulative distribution function of the position of the rightmost point. 
In section|4]we will choose (j) = 0i and (j) = $2, other special cases of pi) , given by 



0i M 



1 for x > 0, 
A for x < 0, 



1 forx>a, 
<h.(x) = {n for0<x<a, 
XjX forx<0, 



to calculate the generating functions ( |2|7[ i at time t 

Wx{x)=H^{x,t), yx^x-a) =H <k (x,t) 

needed to determine the distribution Pn(a) and the density p(a) defined at the end of section]!] 



(34) 



(35) 



3.2 The branching Brownian motion seen from the rightmost point 



The Fisher- KPP equation ( |27| i has two homogeneous solutions: = 1, which is unstable, and = 0, 
which is stable. When the initial condition <p(x) is given by the step function 6(x), see (]32|, the solution 
Hg(x,t) of (27 1 becomes a traveling wave with the phase Hg = invading the phase Hg = 1 II16II21I I81. 



As the front is an extended object, one can define its position m t in several ways; for example one could 
define m t as the solution of Hg(m t ,t) = a for some < a < 1. Here it will be convenient to use the 
following definition 

m t = I dxxd x Hg(x,t). (36) 



One can see using ( 33 1 that m t defined by ( 36 1 is the average position of the rightmost point. 



If the initial condition ( 28 1 is not a step function but is such that (j) (x) = 1 for all large enough x 



and <j)(x) is a constant smaller than 1 for all large negative x, as in (|31|34|, the solution H ( j,(x,t) of (27i 
becomes also a traveling wave. Its position can be defined as in \36\ by 



= I dxxd x H^(x,t). 



(37) 



We are now going to show that the whole measure seen from the rightmost point can be written in terms 
of this position m^: one can rewrite |lT[ i as 



n, x k (zi , ■ ■ ■ ,Zk) = J d* (d x -d Zi d Zk ) y/o A] t ...fe(x, x +zi , . . . ,x +z k ), 

= Yo,l u -M x ' x+Zl '---' x+Zk ) x __ 
- j Ax (a, H \-d Zk ) \i/ QM _ : i k (x,x + zu...,x + z k ), 

= l + { d ;i^ ^ d z k ) / dxxd x y 0i x lr ,.,) it {x,x+zi,...,x + Zk) 



(38) 
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Then from ([30) and ( 37 1 one gets 

%X X X k {zi , ■ ■ ■ ,Zk) = 1 + (d zi H h d Zk 

where is the function pi) with Aq = 0. 



,(*) 



(39) 



Therefore, with the definition ( p7| ) q/f/ie position of the front, the whole information about the measure 
in the frame of the rightmost point, at any time t , can be extracted from the dependence of ntf . 



3.3 The limiting measure and the delays 



In the long time limit, it is known [8.31 1 that the traveling wave solution Hg(x,t) of (27 1, with the initial 
condition (32 1, takes an asymptotic shape, F(x). This means that 

H e (m t +x,t)-^-±F{x), (40) 

where F(x) satisfies 

F"+2F'+F 2 -F = 0, F(-«)=0, F(°°) = 1, J 6xxd x F[x) = 0. (41) 

It is also known, since the work of Bramson [8 1, that, in the long time limit, the traveling wave moves at 
a velocity 2 and that its position ( 36 1 is given by 



m t = 2t — — hit + Constant + o ( 1 ) 



(42) 



If the function fj>(x) is not the step function but is of the form (31 34 1, the solution H<p(x,t) of (27 1 

becomes also a traveling wave with the same shape F(x). This wave is centered around the position w, , 
defined in ( |37j ), and one has 

H^m^ +x,t) — *F{x). (43) 
For large times mf^ is still given by (42 1, but with a different constant |8|. This means that 

m f -m f W >f[<j>], 



(44) 

where f[<j>] is the long time delay in the position of the front due to the modified initial condition, as 
compared to a front starting with a step function. Taken together, ( |43] > and ( p4| ) give 

H^(m t +x,t) >F(x + f[<j)]). (45) 



Using ( |30| >, this becomes 

Vl a ,....,A, k ( m t+x,m t +x + zi,...,m t +x + Zk) 



■+F(x+f[4]), 



(46) 



which shows that the measure of {X\ (t) — m t ,X2(t ) — m r , . ..} (the rightmost points in the branching Brow- 
nian motion seen from the m, frame) does converge when t — > °o to a limiting point measure characterized 
by the functions F(x) and f[(j>]. 

The measure of {^(0 — Xi(t),Xi(t) — Xi(t), . . .} (the rightmost points in the branching Brownian 
motion seen from the X\ (f ) frame) also has a well-defined limit when t — > °°. Indeed, using ( 39 1 and ( |44| , 
one gets 



Xki x k (zi , • • • > Zfc) - — > 1 - (<?-, + • ■ • + ^ J /M 



(47) 



Therefore, in the long time limit, all the information on the distribution of the rightmost points seen from 
X\ (t) is contained in the (j) dependence of the delay f[(j>]. 

Note that, in contrast to ( 39 1 which requires the position to be defined by ( 37 i, the delay f[(j>] in ( |46*1 ) 
or (47 1 depends only on (j): it would not change if we had chosen another definition of the front position. 
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3.4 The superposability of branching Brownian motions 

Let us now consider M independent branching Brownian motions starting at t = at positions u\ , . . . , um- 
Following the same argument as in section [3~Tj the generating function |9]) of the union of the points at 
time t of these M branching Brownian motions is given by the following generalization of ([30| 



M 

WXo,...,X k (x>X+zu---,X+Zk) = Y[ H $( X - U a,t), (48) 



where H&{x,t) is the same solution of (27i with the initial condition (31 1 as in the case of a single 



branching Brownian motion starting at the origin. In the long time limit, using ( |45) , 

M 

VAo \ t {m, +x,m t +x + zu ■ ■ ■ ,m t +x + z k ) F (* + /[#] - u a ). (49) 

r ~ s "* > a=l 

This means that here again, there is a limiting measure when t — > °° for the rightmost points in the m t 
frame. This measure is not the same as before (when one starts with a single point at the origin), as can 
be seen by comparing (j49j and ( |46| ). In particular, the distribution of the rightmost point is different. 
In the frame of the rightmost point, however, one can see using d49l and p8j) that 



%M h (zi . • • • 1 - (^i + • • • + d z k ) fft] > (50) 

as in ( |4"7) >. 

It is remarkable that the generating function % does depend neither on the number M of starting 
points nor on their positions u a . The picture which emerges is that if we superpose the rightmost points 
of several branching Brownian motions, starting at arbitrary positions, the limiting measure in the frame 
of the rightmost point is, when t —¥ °o, the same as for a single branching Brownian motion. 

We will say that, in the long time limit, the measure of the distances between the rightmost points 
in a branching Brownian motion becomes superposable: the union of two (or more) realizations of the 
process (even moved by arbitrary translations u a ) leads to the same measure in the frame of the rightmost 
point as for a single branching Brownian motion. 

As a remark, it is easy to check, that the Poisson process with an exponential density r(x) = e~ ax , see 
section [2.3.2| is an example of a superposable measure: the superposition of M such Poisson processes 
translated by arbitrary amounts u\, . . . , um is identical to a single Poisson process with an exponential 
distribution translated by a -1 ln(e a "' H \-e a " M ). One can also check that, for the same reason, all the 



decorated measures of section 2.3.3 are superposable when Vi is a Poisson process with an exponential 
density. 

In section[6| we will state a stronger version of the superposability property of the branching Brownian 
motion. 



4 Some quantitative properties of the branching Brownian motion seen from the rightmost point 



In this section we obtain, by integrating numerically the equation ( 27 1 with the appropriate initial condi- 
tion, some statistical properties of the limiting measure seen from the rightmost point. 



1 1 



4. 1 Average distances between consecutive points 



The analytic calculation of the delay f[<j)] is in general not easy. For <j> = <j)\ given by (34i, however, it 
was possible to show [9| that when 1 — X <C 1, the delay is given by 

/[ft] ~ -ln(l-A) -ln[— ln(l-A)] + 0(1), (51) 

and, from this, one could deduce that, in the long time limit, the average of the distance between 
the z'-th and the (/+ l)-th rightmost points is given for large i by 

(d u+l }~--^-.. (52) 
i i mi 

In |9j, the numerical values of the distances between the rightmost points were also obtained by inte- 
grating the Fisher-KPP equation with the initial condition <j>\ in ( 34 1 and by using (|6]i (in practice we 



integrated numerically the equations satisfied by the coefficients of the expansion of y/^ (x) in powers of 
X). It was found that 

(<* 1>2 )~ 0.496, (</ 2 , 3 ) -0.303, (rf 3 , 4 ) ^ 0.219, 

(d 4 . 5 ) -0.172, (</ 5 , 6 ) -0.142, (d 6J ) -0.121. <Xl) 
The results ( |52|53[ > gave evidence that the distances between the rightmost points of the branching Brow- 



nian motion were different from those of a Poisson process with an exponential density (21 



4.2 Distribution of the distance between the two rightmost points 



According to ( 13 1, to obtain the distribution Paifl) of the distance between the two rightmost points, one 
needs to calculate y/o^(x,x — a) to first order in jj,. We first remark that at time t for a > 0, \\foa(x,x — a) = 
Proba[n(x) = and n(x — a) = 0] = Proba[n(x — a) = 0] = Hq(x — a,t) where He(x,t) is the standard 
Fisher-KPP front with the step initial condition (it is also easy to see from the definition ([34} of 02-) Then 
writing at time t that 

Yon[x,x — a) = Hg(x - a,t) + jj.R a (x — a,t) + ffijJ 1 ) (54) 
is solution of the Fisher-KPP equation, and using the initial condition 02 m 1 )34) , one gets 

d t R a = d?R a -R a +2H e R a ; R a (x,0) = \ \ fo ^ fl ^ x<0 ' (55) 

I otherwise. 

Then, from ( fT3j ) one gets 

P 12 (a) = -d„ 2 fdxR a (x,t). (56) 



Figure [3] shows our numerical result for the distribution Pn(a) of the distance between the two right- 
most points in the long time limit. More details on our numerical procedure is given in appendix |A| 

We see that P\2{a) is very close to 2e~ 2a for the values of a which have a significant probability of 



occurring. This is of course consistent with an average distance (53 1 close to 1 /2. For large a (events with 
a small probability), however, the exponential decay is faster. We now present a simple argument leading 
to the following prediction, which is consistent with our numerical data, 

P n (a) ~ e - (l+ ^ )a for large a. (57) 

In the long time limit, the right frontier of the branching Brownian motion moves at velocity v = 2. Let us 
assume that a large distance a between the two rightmost points is produced by the following scenario: by 
a rare fluctuation, the rightmost point escapes and, without branching, goes significantly ahead while the 
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0.5 1 1.5 2 2.5 3 50 100 150 200 250 300 350 400 

a a 

Fig. 3 The density of probability P\i[a) of observing a distance a between the two rightmost points in the r — > °° limit, as a function 
of a. For a small (left part), the distribution is very close to 2e~ 2a . For larger values of a, one observes a faster exponential decay of 
order e -(l+V5)a 



rest of the points go on as usual, forming a frontier moving at velocity v = 2. Such an event leads to the 
distance a between the two rightmost points if, during a time T, the rightmost point moves (by diffusion 
alone) by a distance a + 2t without branching. The probability of such a scenario is 



Proba(Xi — X2 ~ a after an escape time t) ~ exp 



(a + 2t) 2 
4t 



(58) 



The first term is the probability of diffusing over a distance a + 2t during time T, and the second term is 
the probability of not branching. The probability to observe a large distance a 3> 1 is then dominated by 
the events with % chosen to maximize d58ll, that is 



^optimal 



2V2' 



(59) 



and this leads to (57 1 in good agreement with the numerical data of figure [3] 

There is a remarkable relation between the decay rate in ( |5"7] > and the shape of the traveling wave 
solution of (41 1. Around the stable region F = 0, the equation (|41|i can be linearized and one has 



F(z) ~ Cste x e n for z ■ 



with 



■2r-l =0. 



(60) 



We emphasize that this is a linear analysis of the stable region, which is usually uninteresting (in contrast 
to the unstable region which determines the velocity). The solutions for r are 



-1±V2. 



(61) 



r = — 1 + \/2 is the correct root as it is the only positive solution and F(— °°) has to vanish. The other 
solution r = — 1 — v2 (the wrong root) coincides (up to the sign) with the decay rate of the distribution 



Pn{a) for the distance a between the two rightmost points ( 57 1 



As explained in appendix [Bj this coincidence exists in a broad class of branching processes: each 
variant of the branching Brownian motion is linked to a variant of the Fisher-KPP equation, and the 
wrong root in the linear analysis of the stable region always gives the asymptotic decay rate of P\2{a). 



13 



4.3 Average density seen from the rightmost point 
To obtain the average density of points at a distance a on the left of the rightmost point, one needs, 



according to ( |14) , to calculate \^q^{x,x— a) for jj. close to 1. As in section 4.2 we first remark, from the 
definition |F7), that y/oi (x,x — a) = Proba[n(x) = 0] = Hg(x,t) is the standard Fisher- KPP front with the 
step initial condition. Then, writing at time t that 

Yon(x,x-a) =H g (x,t) - (l-n)R a (x,t) + - fi) 2 } (62) 

is solution of the Fisher-KPP equation, and using the initial condition 02 in < [34] >, one gets 

d t R a = d 2 R a -R a +2H e R a - R a {x,0) = \ 1 for0 ^ x<a > (63) 

I otherwise. 

It is the same equation as for R a in ( [55) , but with a different initial condition. Then, from ( fl"4| ) one gets 

p{a)=dj JdxR a (x,t). (64) 



a 



One can notice the great similarity between the expressions for the average density p(a) of points at 
distance a from the rightmost ( 63|64 i and the probability distribution Pn(a) for the distance between 



the two rightmost points ( |55|56[ ): one goes from one to the other by simple changes of signs, as in the 
example of a Poisson process (|18|. 

Figure |4]presents our numerical results for p{a) in the long time limit. We see that p(a) increases as 

p(a) ~ Cste x ae a for large a. (65) 

Note that a Poisson process with such a density would lead to asymptotic distances between points given 
by ( |52) . The branching Brownian motion is however not a Poisson process as the points are correlated, at 
least near the tip. 




Fig. 4 The average density p(a) of points at a distance a of the rightmost in the long time limit grows like ae" . When the data is 
multiplied by e~", as shown in the figure, the linear prefactor is clearly visible. 
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5 Generalizations to other branching processes 

All the results of sections [3] and |4] can be generalized to other branching processes on the line where 
points move and branch independently of the positions and of the motions of the other points. In such 
systems, the function H^(x,t) defined in ( |25| l is also solution of an equation similar to the Fisher- KPP 
equation ( p7| ). Here are four examples: 

A) The points perform Brownian motions and branch as before, but at each branching event there is a 
probability p to branch into three points and 1 — p to branch into two. Then H,j,(x,t) evolves according 
to 

d,H$ = 9 2 Ha +pHl + {\- p)B\ - Ha . (66) 

B) Time is discrete with steps of duration 8; at each time step, a point at position x branches into two 
points at positions x + £\ and x + £2, where the e, take independent random values distributed according 
to some given p(e). The evolution of HA(x,t) is then given by 



dep(e)H^(x-e,t) 



(67) 



In this example, the positions of the points can be thought of as the possible energies of a directed 
polymer on a Caley tree with independent random energies e on the edges of the tree B 1411271 . 

C) Time is continuous but space is discrete with steps 1; during dt, each point at position x has a prob- 
ability df of being removed and replaced by two points at position x+l. The equation satisfied by 

H$(x,t) is 

d t Hf (x, t)=H^x-l,t) 2 -H^(x,t). (68) 

This example is relevant to the theory of binary search trees 128 . 2411231 . 

D) Time and space are discrete with steps s for space and 8 for time; in a given time step, a point at 
position x has a probability 8 of branching into two points at position x, a probability 8/s 2 of jumping 
to the left, 8/s 2 of jumping to the right, and 1 — 5 — 28/s 2 of doing nothing. Then: 



fy(x,t + 8) =H (j) (x,t) + 8 



Ha (x - s, t) +Ha (x +s, t) - 2Ha (x, t) 



-H^{x,t)+H^{x,t) 



(69) 



which is of course a discretized version of the original Fisher-KPP equation. ( |69| > is actually the 
equation we used in our numerical simulations, see appendix [A] 

In all cases, these equations have Ha = 1 as an unstable fixed point, and Ha = as a stable fixed point. 
For initial conditions Ha(x,Q) = <j>(x) of the type (31 32 34 1, the function H<j,(x,t) develops into a trav- 
eling wave moving at a specific velocity v* . We recall briefly the procedure to determine the asymptotic 
velocity v* of the front (which is also, through ([33|, the velocity of the rightmost point in the branching 
process). One looks for traveling wave solutions moving at velocity v of the form HA.(x,t) — F(x — vt) 
and solve the linearized equation around the unstable fixed point by writing 1 — F(x) ~ Ee~Y x . This leads 
to a relation between 7 and v, and the minimal value v* of v reached at some 7* is the velocity selected 
by the front [3 1 ']. (We only consider here cases where the function v(y) has a minimum.) For our four 
examples 



A) v = 7+^:v< 



C) 



y 

lel-\ 



:4.311, 



P, 



B) v=^ln 
7<5 

D) v^-^ln 
70 



2 J Aep{e)e y£ 



^cosh(7*)-l +g 



(70) 



Once the equation for Ha of a particular branching process is written, one has access to all the gen- 
erating functions i//)^(x,jic — a), etc., see ( 2 7|9 1, by choosing the appropriate initial conditions 
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( 31|34b for the front equation. The whole measure in the frame of the rightmost point is then obtained 



from (|39|l at any finite time t . Note that to prove the existence of a long time limit to the point measure 



in this frame for a specific branching process, one would need a version of Bramson's result ( |42| i for this 
process which is, to our knowledge, not known in the general case. 

It is natural to ask which properties of the branching Brownian motion can be extended to other 
branching processes. If the measure for the distances between the rightmost points has a long time limit, 
then the arguments of section [3~4] can be easily generalized and one can show that it is superposable. We 
have checked that the analytical argument |9| leading to the asymptotic expression |52]i for the average 
distances (rfjy+i) at large times can be extended in case B) for a large class of densities p(e) and yields 

(di,i+ 1 ) — ~~ ( ~ Tj~ '. ] for large i. (71) 
Y \i i In i J 

We have also checked numerically on examples C) and D) that the density at a distance a of the rightmost 
point is, as in ( [65) , 

p (a) ~ Cste x ae 7 *" for large a. (72) 

For the tail of the distribution Pn{a) of the distance a between the two rightmost points, we discussed a 
scenario, at the end of section[42j which can be generalized (see appendix[B]) to calculate the exponential 
decay Pn(a) for more general branching processes. This scenario, however, can only hold if points can 
move without branching, as in our examples A) and D); for instance, in example A), it predicts an expo- 
nential with a decay rate equal to ^/l + p + y/2 + p. In examples B) and C), the points branch whenever 
they move and the tail of Pn(a) is in general not an exponential. 

Note that special care should be taken if the points are located on a discrete lattice, as in cases C), D) 



and possibly B): quantities such as P\2{a), see ( 13 1, become probabilities rather than densities of proba- 
bility and quantities such as p(a), see ( 14 1, become average numbers rather than average densities, and 
all the formulas in the previous sections need to be adapted: integrals become discrete sums, derivatives 
become finite differences, etc. If one interprets n(x) as the number of points strictly on the right of x, then 



the generating functions l/fy, y/X^, . . . are still related to as i n ([30 35 1 with the choices ( 31|34 1 for the 



initial condition <j). Then, for instance, one can show easily that (|6]i becomes 

Y t x[Yo{x)-Yo(x-s)-^x{x) + Yx(x-s)] = £{d i j + i)X i , (73) 

x i> 1 

where s is the lattice spacing. For these systems on the lattice, there are new properties that can be 
investigated. As an example, if N is the number of points on the rightmost occupied site, then it is easy to 
check that 

</*"> = 1+1 [y% (x,x - s) - H e (x, t)] . (74) 

X 

The whole distribution of N can then be determined by numerical integration. In the case of our ex- 
ample C), the number N corresponds to the number of leaves at the deepest level in a binary search 
tree ll28ll and we found numerically that, at large times, Proba(A r = 2) ~ 0.50, Proba(A r = 4) ~ 0.23, 
Proba(A? = 6) ~ 0.11, Proba(N = 8) ~ 0.06, etc. 



6 Large time measure in the frame of Lalley and Sellke 



We recall from the results of section [371] see ( |40] i and ( p3) l, that the distribution of the rightmost point is 
given in the long time limit by 

limProba[X"i(?) <m t +x] =F(x), (75) 
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where m t is the average position of the rightmost point and is asymptotically given by ( |42| i and where F 
is defined in ( |4Tj ). The property ( |75) l is however non-ergodic, as shown by Lalley and Sellke [22 1 in the 
sense that for a single realization of the branching Brownian motion 

[fraction of the time that X\(t) <m t +x] ^F(x). (76) 

In fact, the left hand side of f76] > is not self-averaging and depends on the realization. This is illustrated in 
figure [5] for the realization on the left, there were few branching events at early times and the first points 
wandered to the left, leading at larger times to an asymmetric picture. For the realization on the right of 
figure [5] there were many branching events early. For the right realization, the rightmost point is almost 
always on the right of m t while it is almost always on the left of m t in the left realization. 




Fig. 5 Two realizations of the branching Brownian motion up to time ( = 20. The horizontal direction represents space, and time 
increases downwards. The dotted gray lines are m, and — m t , the average positions of the rightmost and leftmost points, as measured 
from equation j36fr . 

Visually, these strong memory effects of the early stages of the branching Brownian motion do not 
seem to decay with time, and it looks like the fluctuating right frontier of the system settles at some 
random fixed distance C from m t : 

X l (t)=m t +C + T]i(t) for large f, (77) 

where C would depend on the realization but not the time, and where T]i(t) would be a time-dependent 
random number centered around zero. A natural question is whether it is possible to define C for each 
realization in such a way that the distribution of ?7i(f) becomes in the long time limit independent of C 
and t, the idea being that the branching Brownian motion at long times seen from m t + C would "look the 
same" for any realization, whatever is the value of C. 

A related question was addressed by Lalley and Sellke [22] in the following way: for each realization 
of the branching Brownian motion, define Z as 

Z = limZ, where Z, = £[2f -X i (f)]e* i W -a (78) 

i 

(The sum is over all the points X{(t) present at time t .) As shown in 1221 . Z t has a limit Z for almost every 
realization; that limit is finite and positive. Lalley and Sellke prove then a limit theorem for the frontier 
of the branching Brownian motion which we interpret as follows: 

limProba(X,(?) <m,+x\z) = exp ( - AZe'*) = exp ( - e"^ ln(AZ) l) , (79) 
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where A is a constant related to the large x behavior of F(x), see ( |101| i. In words, if one considers only the 
realizations of the branching Brownian motion with a given value of Z, then the large time distribution of 
the rightmost point is given by a Gumble located around m t + ln(AZ). To make the link with ( |77| , 

C = ln(AZ), Proba(77i < x) = exp ( - e~ x ) . (80) 

In appendix[Cj we present the precise theorem stated by Lalley and Sellke and we argue that |79) should 
be equivalent to their result. 

A natural extension to Lalley's and Sellke's result is to write for all the points ; 

Xi(t)=m t +]n(AZ) + r]i(t), (81) 

as in ( |77| ), and ask whether the joint distribution of T[i(t), T\2{t), T73 (?), . . . reaches a long time limit which 
is independent of the value of Z. We show in appendix [C] that Lalley's and Sellke's result can be extended 
to all the generating functions . Our interpretation of this extension is 

lim(n<K OT <+ x -^( f ))| Z ) =exp(-AZ e -- v - / W) = exp ( - e^- 1 "^^]) , (82) 

i 

where the delay function f[<j>] is the same as in ( |45] >. By choosing = (the step function), ( |82"| ) reduces 
to |79| ). By choosing as in ( fJTj ), one sees from ( [82) that the distribution of points at the right of the 
branching Brownian motion conditioned by Z reaches a long time limit where Z only appears through the 
global shift ln(AZ). This means that at large times, the distribution of the rightmost points in a branching 
Brownian motion has a well defined measure independent ofZ located around m, + ln(AZ). 

As an example, if one chooses the function (j>\ defined by p4| ), one can easily show from ( [82") ) and ( |5Tj ) 
that, in the m t + ln(AZ) frame, the average density of points at any position is infinite in the long time 
limit. 



6. 1 Superposability property 

If one considers two branching Brownian motions a and b starting at arbitrary positions, then the points 
in a at large time will be characterized by a random value ZW and a realization of the point measure 
described by ( f82] >; idem for the points in b. If one considers the union of these two branching Brownian 
motions, one gets from (|82|) 



lim(jl0(jn, +*—$(»)) ZW,Z^) =exp(-AZe-- v - / W) =exp(-e-[ x - ln ( AZ ) +/ W]), (83) 



with Z = +zW. This means that the point measure reached in the long time limit in the m, +ln(AZ) 
frame is the same whether one started initially with one, two or, by extension, any finite number of initial 
points at arbitrary positions on the line. What does depend on the initial number of points is only the law 
of the random number Z, not the positions around m t +ln(AZ). This is to be related to the discussion in 
section 3.4 where we showed that, in the long time limit, the measure seen from m, depends on the initial 
number of points while the measure seen from X\ (f ) does not. 

Furthermore, the large time measure of the points in the m t +ln(AZ) frame has the following property: 

Starting with two realizations } and {T},- }, then for any pair of 

real numbers a and /3, the ensemble of points {r\\ a "' + a} U {tj ; . + j3} (84) 

is another realization of the same measure shifted by ln(e a +e&). 

(Think of {tj^ } as the offspring of a in the m t +ln[AZ' a ^] frame and of a as InZW; idem for b. The 
shifts a and /3 are arbitrary because Z" and are unbounded independent random numbers.) 
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The property ([84} of the point measure in the m, + ln(AZ) frame is a stronger version of the super- 
posability property discussed in section [34] clearly, it implies that the distribution of distances between 
the rightmost points is invariant by superposition, but it gives more information on the measure as it en- 
compasses the position of the rightmost point. In particular, one can check that, in any such measure, the 
rightmost point is Gumble distributed. 

The simplest point measure with the property (|84"|) is the Poisson process with an exponential density 



Ke~ x , for an arbitrary K. Furthermore, all the decorated measures of section 2.3.3 when Vi is a Poisson 
process with an exponential density Ke~ x are also superposable measures. A natural question is then: 
can any superposable point measure be constructed as a decorated exponential Poisson process for a 
well chosen decoration measure ? A more specific question would be: is the limiting point measure of 
the branching Brownian motion in the m t +ln(AZ) frame can be constructed as a decorated exponential 
Poisson process? 

We conjecture that the answer to the last question is yes. Indeed, it is known |[T4l [3l that if one consid- 
ers two points chosen at random among the rightmost points in a branching Brownian motion at a large 
time t , then the time one needs to go back to find the most recent common ancestor of those two points 
is either very short (of order 1) or very long (of the order of the age t of the system). This means that one 
can group the rightmost points into families where two points belong to the same family if the branching 
event that generated them occurred recently. The branching Brownian motion could then be a decorated 
exponential Poisson process, where the rightmost points of all the families are Poisson distributed with 
an exponential density |3| and where the members of a given family form the decoration. This interpre- 
tation helps to understand a question raised by our previous work |9|: the distances between points in the 
branching Brownian motion are given by ( 52|53 I, but if one considers only the rightmost point of each 
family, then the distances are given by \2\\ with a = 1, as in the GREM or in the Ruelle cascade, see 
section l2.3.4l 



7 Conclusion 

In the present work, we have shown that all the statistical properties of the rightmost points in a branching 
Brownian motion can be obtained by solving a front equation with a well-chosen initial condition. The 
distribution of the positions of the rightmost points seen in the frame m t (the average position of the right- 
most) has a long time limit. The properties of the limiting distribution can be expressed as the long time 
delays of the traveling wave solution of the Fisher-KPP equation when one varies the initial condition. 
This limiting distribution is however modified if one considers the union of several branching Brownian 
motions. 

If one considers, however, only the distances between the points, for example if one looks at the 
distribution of all the positions of the rightmost points seen in the frame of the rightmost one, one obtains 
at large times another limiting distribution which does not depend on the initial positions of the branching 
Brownian motions (as long as there are finitely many of them). We called this property superposability. 

In section [4] we have measured a few properties of this limiting distribution, and in section [5] we 
explained how our results can be extended to more general branching random walks. Lastly, in section[6] 
we argued that in the Lalley and Sellke frame, the branching Brownian motion satisfies a stronger version 
of the superposability property, see (84) . 

In the future, it would be interesting to characterize more precisely the limiting measure of the 
branching Brownian motion and of the branching random walks to see whether some universal prop- 
erties emerge. For example, we believe that the average density seen from the rightmost point should 
always grow as in f72"] i. It would also be interesting to check whether, in Lalley's and Sellke's frame, the 
limiting distribution is indeed a decorated exponential, and to determine the properties of the decorating 
measure. 
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The question of ergodicity, raised at the beginning of section|6] is also an interesting open question. 
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A Numerical simulations 

We performed the numerical simulations by discretizing the branching Brownian motion in space (with a grid length 8) and time 
(with steps i) as in the example D) of section [5] The corresponding front equation is given by {69) , which is of course the most 
straightforward discretization of the Fisher-KPP equation (27) . The solutions to the discrete equation (69) converge to the solutions 
of the Fisher-KPP equation (27) , if s — > with 8/s 2 held constant and small enough. 

We used three sets of values for s and 5 and computed the exact asymptotic ve loci ty v* of the front and the decay rate y* 
of the asymptotic shape F[x) by minimizing numerically the function v(y) given in (70p ). The decay rate j8 for the probability 
of observing a distance a between the two rightmost points, see equation {57}, was computed using the recipe given at the end of 
section |4~2] and explained in appendix|B] All these values are presented in table[T| 



s 8 


v* f P 


0.25 i 2 /5 = 0.0125 

0.10 .v 2 /5 = 0.002 

0.05 .v 2 /6 ~ 0.000416667 


1.980480133 1.004581693 2.387337826 
1.996840367 1.000747277 2.409772891 
1.999375296 1.000104046 2.412897517 


Fisher-KPP 


2 1 ^2+1 =2.414213562 



Table 1 Values of v*, 7* and p for our discretized branching processes, compared to the values in the Fisher-KPP case. 



The simulations were made on a finite but large domain centered around the position of the front; typically it extended to about 
a distance 1000 ahead and behind the center of the front (respectively 4000, 10000 or 20000 lattice sites depending on s). The 
values at x = ±°° were exactly computed and used for the boundaries of the domain. Whenever the front moved by more than one 
unit space lattice, the whole data set was recentered. The simulations were performed up to large times of order 10000 (0.8 to 24 
millions of time steps) and the data was extrapolated to obtain a value at t = °°. To do this extrapolation, we used a more precise 
asymptotic expansion of the position of the front than j42j: according to 1151 . 

, 3 Q1/2 at a 3/2 
m, = v*t- — In t + Constant + — ^ H h -Ar + --- , (85) 

where the number does not depend on the initial condition. As we measure the delay m, — m, , many terms cancel and one 
gets 

m ,- m ^ =m+ S ^ + S ^ + .... (86) 

All the quantities we measure are derivatives of see (47) , and have therefore the same large time expansion as (86) . Thus, we 
extrapolated our numerical data to the large time limit by fitting it with the function A +B /t + C/t 3 l 2 for times larger than (typically) 
5000, see figure|6] and by using A as the end result. 

On figure [3] the three data points were presented together; on figure [4] we have drawn together for each data set the function 
p(Xi — a)e~"f a using in each case the value of 7* of table[l] In both cases, the superposition was nearly perfect, and so we expect 
that on the scales of the figure, the curves would not change noticeably for smaller values of s and 5. 



B Distribution of the distance between the two rightmost points 

In this appendix we generalize, to any branching random walk, the argument leading to the asymptotic decay (57) of the distribution 
of the distance between the two rightmost points in the branching Brownian motion. 

We consider a generic branching random walk in discrete space (with spacing s) and time (with intervals 8) defined by the 
following family of functions 

. \ — ( T ne probability that a point at position x branches during a time 

Pn\ 1 1 ■ ■ ■ i r n) y s i e p i n to n points located at positions X + r\, ...,x + r n . 
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Fig. 6 Average density at a distance 30 from the rightmost point as a function of time (symbols), fitted for ; > 5000 by the function 
A + B/t + C/t 3 / 2 (line). The inset shows the quality of the fit by displaying the ratio between the data points and the fitting function. 



We assume that po = 0, so that there is no extinction. Then p\ (r) can be thought as the probability that the point does not branch 
but moves by a distance r. The continuous time and/or space cases can be obtained as suitable s — > and/or 5 — > limits. 

Let exp[fg(j8)] be the generating function of the displacement during time 1 of one point conditioned on the fact that this point 
does not branch: 

gW(p) = ^e' 5 ' Proba(the point moves a distance r without branching during time t). (88) 
As the time steps are independent, the function g(p) can be computed during the time interval 5 which gives 

*»«<» (89) 

r 

Note that g(0) < as soon as the branching probability is non-zero. We want now to evaluate the probability that a point moves a 
distance r, without branching, during time t. For large t, it takes the form 

Proba(the point moves a distance r without branching during time t) ~ exp \tf y- J J , (90) 

where f(c) is a large deviation function. Using (88) , one finds that /(c) and g(j3) are related by a Legendre transform 

g(j3)=/(c)+j3c. 

Now, assuming as in section |4~2] that the events which contribute most to a large distance a between the two rightmost points 
are those where the rightmost point moves, without branching, a distance a ahead of the frontier of the branching Brownian motion, 
one gets 



P\2{a) ~ max < exp 



z 

where v* is the velocity of the front. For large a, the optimal T is also large and it satisfies, by derivation, 



(92) 



Let c = a/z + v*. Using (91) , equation j93} becomes 

g(j3) = j3v*. (94) 
Remarkably, this equation does not depend on a. Replacing into j92j gives 

P n (a)~e-V a . (95) 



The asymptotic decay rate of the probability distribution function of the distance between the two rightmost points is therefore 
simply the positive solution j8 of (94) with g(j3) given by j89}. In the branching Brownian motion, g(j8) = /3 2 — 1, v* = 2, so that 
(94) gives indeed j8 = 1 + \/2. 

As can be checked easily from {89| , the function g(/3 ) is convex. Therefore, as g(0) < 0, equation (94) has at most one positive 
solution and at most one negative solution. The positive solution is the relevant one here. 
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We are now going to show that the negative solution of (94) , if it exists, gives the asymptotic shape of the associated traveling 
wave. We write the front equation associated to the branching point process (87) . Using the same method as in section[3] we find 

H t (x,t + S) = '£ E Pn(ru---,rn)flH^x-ri). (96) 

n>lr\ r ..,r n i=l 

We look at the shape F of the traveling wave solution for a step initial condition, which moves asymptotically at the velocity v* : 

H^{x,t)=F{x-v*t). (97) 

Using (96) , we see that in the stable region F(x) -C 1 one has 

F(x-v*S) =Y d Pl(r)F(x - r) + ff(F 2 ). (98) 
r 

We look for an exponential solution to this linearized equation: F(x) ~ e** with A > 0, as F(— °°) = 0. (Note that a periodic 
modulation of this exponential could occur as r takes only discrete values.) Inserting into (98) and using (89) , one finds that the 
equation for A is 

g(-k) = -Xv*, (99) 

which is the same equation as (94) for A = — /3 . 

To summarize, a positive solution to (94) gives the exponential decay rate of the probability distribution of the distance between 
the two rightmost points, see (95) , while a negative solution gives the coefficient —A governing the shape of the front F{x) in the 
stable region F(x) -C 1. 



C Lalley's and Sellke's result 

Lalley's and Sellke's theorem |22| is 

lim lim Proba (X, (t) < m, +x\ {Xds)}) = exp ( - AZe~- r ) , (100) 

S — leaf — ► V 

where A is the constant appearing [8] in the large x expansion of the function F(x) defined in (41) 

F(x)~l-(Ax+B)e- x for large*, (101) 

and Z is defined in (78) . In words, given the positions {X, (.v)} at time s, there is a t — > °° limit to the probability that the rightmost is 
on the left of m, +x which depends, obviously, on the {X, (.s)} and is as such a random variable. As s goes to infinity, this random 
variable converges almost surely to the Gumble distribution around ln(AZ) . 

This result can be extended into the following: for any suitable function (see section |3~2) , one has 

lim lim (T[$ [m, +x-Xi(t))\{X t (s)}\ = exp (-AZe-^^) , (102) 
f' 

where f[(j>] is the delay function (45) . For (j> = 8, (102) reduces to (100) . 

We first give an outline of Lalley's and Sellke's proof applied to the case (102) . Given the positions Xj (.v) of the points at time s, 
the system as time t > s can be seen as a collection of independent branching Brownian motions at time t — s starting from the Xj(s). 
Therefore 

(n*[*-X.-(0] {*>(*)}) = YlH4,(x-Xi( S ),t-s), (103) 
( i 

where the product in the right hand side is made on all the points present at time s. 

We replace x by m, +x, to center around the position of the front, and suppose ; large. It is easy to see from Bramson's formula 
(42) that m, = m,- s + 2s + o(l) as t becomes large, so that 

(jltim, +x-X,(t)] I = (m,- s + 2s + x-X i (s)+o(l),t-s) , (104) 

i i 

and, using (45) , 

lim (]J<l>[m t +x-X,(t)] I {*,(*)}) = Y[F(2s+x-X,(s) +/[(/>]) . (105) 

i i 

We now take s large. Of all the points present at time s, the rightmost is around 2.v— | logs, see (42) . Therefore, 2s — Xi(s) diverges 
for all i. Using (TTJI) , 

lim (J}* fa +*-*(')] |{*/W>) ^exp(-£ [A(2s + x- Xds) + f [(!>}) +B] e^-^ 'W-^l) . (106) 
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Following Lalley and Sellke, we introduce the quantities 

Y s = ^e- 2s+x ^, Z s = Y,[2s-M*)]e- 2s+XM , (107) 
f i 

see (78) , so that 

Jm^Y[<p[m, +x-Xi(t)]\{Xi(s)}} ~ exp (- [AZ S + (Ax+Af[(j>] +B)Y S ] e^M) . (108) 

i 

Finally, the most technical part of Lalley's and Sellke's proof is that Y s and Z, are martingales converging when s — > °° to lirn,-^ F s = 
and limj-^Zj = Z > respectively, which lead s to (102) . We do not reproduce this part of the proof here as it does not concern 
our extension with the function <j> and it works in [102) exactly as in (100) . 

In (102) , the average is made on all the realizations with a given set {X, (s)} of points at a large time s but the only relevant 
quantity appearing in the generating function (108) is Z,. One would obviously have reached the same result if one had conditioned 
by Z s instead of by the {X, (.v)}. Furthermore, as Z, converges quickly to Z, as illustrated on figure[5] we argue that conditioning by 
Zj at a large time s or directly conditioning by Z should be equivalent, hence (79|82) . 
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